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CNl ' Abstract 

>: 

OO . We use our former results on tt + it 5-wave obtained in a nearly 

^ \ assumption-free way from the 17.2 GeV/c data to predict the 7r°7r° S- 

• wave. The predictions are compared with the recent results of the E852 

experiment at 18.3 GeV/c. A good agreement is found for only one (the 
"down-flat" ) solution while the second one (the "up-flat" ) is excluded 
by the 7r°7r° data. Thus the long-standing "up-down" ambiguity has 
been finally resolved in favour of the S-wave intensity which stays large 
and nearly constant up to the KK threshold. A joint analysis of both 
sets of data leads to a reduction of errors for this solution. 
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^ 1 Introduction 

Despite of several decades of efforts a spectrum and properties of scalar 
mesons are still not well known. This is particularly true in the low mass region 
(below the KK threshold), where the whole information comes from the tttt 
partial wave analysis (PWA). It is interesting to note that the best experimental 
data on the 7r + 7r~ system were collected in the seventies. We mean here e.g. 
the CERN- Munich experiment [p, which supplied 3xl0 5 events of the reaction 

7T~p — ► -K + r K^n (1) 

at 17.2 GeV/c. It should be stressed that the number of observables provided 
by such experiment is much smaller than the number of real parameters needed 
to describe the partial waves. Consequently, special physical assumptions^ were 
made in most studies [1-4] of the 7r + 7r~ system and are still made in all PWA's 

1 Hereafter called the standard assumptions. 
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of the 7T 7r° system. It has been customary to produce results without even 
mentioning the assumptions essential for their derivation. 

These assumptions are generally connected with ignoring of the role of nu- 
cleon spin. Also the a± exchange amplitude (which is dominantly spin-nonflip) 
was neglected thus assuming the full dominance of the pion exchange. It should 
be stressed that the reality can be more complicated, the nucleon spin cannot 
be a priori ignored and there is a significant pseudovector (the a± ) exchange in 
addition to the dominant pseudoscalar pion exchange. This was demonstrated 
by another important experiment, namely that of the CERN-Cracow-Munich 
(CCM) collaboration || which yielded 1.2 x 10 6 events of the reaction: 

T[~P] — > 7T + 7T~n (2) 

also at 17.2 GeV/c. The presence of the transversely polarized target provides 
additional observables. It was shown that there is a considerable a\ exchange. 
Combination of the results of both experiments enables a quasi-complete PWA 
with much weaker physical assumptions. The analysis is only a quasi-complete 
one because we do not know the phase between two sets of transversity ampli- 
tudes. Nevertheless e.g. the full (containing both pseudoscalar and pseudovec- 
tor exchange as well as both 1 = and 1 = 2 components for S-wave and 
D-waves) partial wave intensities could be determined in a completely model- 
independent way. This was done in another CCM paper ||. While for dipion 
mass m n7T = (1000 1800) MeV a unique solution was obtained, the up-down 
ambiguity reappeared below 1 GeV. 

In our first paper (hereafter called paper I) we "have constructed a bridge" 
between two sets of transversity amplitudes for the S'-wave. This was done by 
first demanding the phases of the leading P-, D- and F-transversity ampli- 
tudes to follow exactly the phases of the Breit-Wigner p(770) , /2(1270) and 
P3(1690) resonant amplitudes in the low, medium and high mass region, respec- 
tively. Then, from the measured phase differences between these waves and 
the S'-wave we calculated the absolute phases of the S-wave transversity am- 
plitudesQ Finally the pseudoscalar and pseudovector exchange amplitudes in 
the S-wave were explicitly determined from the relevant combinations of the 
S-wave transversity amplitudes. In order to calculate the 1 = components 
we have used the 1 = 2 S-wave as measured by the CERN-Munich || col- 
laboration. This separation, first ever made, has been done using much weaker 
assumptions^ than those made in any earlier analysis. 



2 Let us note that such a procedure would be impossible for the low mass 7r°7r° S'-wave 
even if the reaction (3) was measured on a polarized target. There would be no P— wave to 
be used as a handle to determine the absolute phase of the S'-wave . 

3 The main assumption was neglecting of a possible influence of the a\ exchange in J = 2 
S-wave . This assumption was also made in all other studies. 
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The price we paid was a fourfold ambiguity in our pseudoscalar exchange 
S-wave amplitude. In addition to the "up-down" ambiguityQ of the old CCM 
analysis U there were ambiguities resulting from adding or subtracting the 
phase difference since the PWA of the CCM group yielded only the absolute 
value of the phase difference. Thus we have "down-flat" , "down-steep" , "up- 
flat" and "up-steep" solutions. Differences between "flat" and "steep" refer 
mainly to the behaviour of the 5-wave phase shifts below the / (980) . Above 
the /q(980) resonance all the solutions are fairly similar. The main difference 
is that both " steep" solutions contain a relatively narrowf] resonance under the 
p(770) (like the old "up" solution). Such resonance was persistently claimed 
by Svec [[J. On the other hand, both "flat" solutions indicate a /o(500) state 
[ |TC| ] with a width of about 500 MeV. In particular the "down- flat" solution is 
essentially similar to the old solution of the CERN-Munich group [JJ obtained 
with the help of the standard assumptions. 

Having arrived at as many as four solutions we have immediately started a 
thorough investigation of them in a hope of reducing their number. Already 
in paper I we have shown the "down-steep" solution to violate unitarity while 
both "flat" solutions easily passed this test. The "up-steep" solution although 
behaving in a queer way could not be immediately excluded. Our next paper 
[Pll (hereafter called paper II) was devoted to a detailed study of this suspi- 
cious solution. Using the experimental fact that there is no significant inelastic 
channel below the /o(980) we have imposed the r\ = 1 inelasticity for all mass 
points in this mass region. For the "up-steep" solution this led to unphysical 
results in 7 out of 20 mass bins and queer behaviour in the remaining mass bins. 
The effect was even stronger for the "down-steep" solution, already rejected 
in paper I. Thus removing both "steep" solutions we were left with two "flat" 
solutions equally well describing the 7r + 7r~ data. The further step demanded 
an extra input. 

Fortunately such an input has been recently provided by the BNL E852 
collaboration ||12|| , which has collected 8.5 x 10 5 events of the reaction 

n~p — > 7r°7r°n (3) 

at 18.3 GeV/c. A comparison of both "flat" solutions to the data is a 

topic of the present paper. 

The paper is organized as follows. In Sect. 2 we describe the 7r + 7r~ and 
7r°7T data. In Sect. 3 we relate our 7r + 7r~ 5-wave amplitudes to the amplitudes 
describing the 7r°7r° channel, in which the role of the isospin 2 component is 
substantially enhanced in comparison with the 7r + 7r~ channel. In Sect. 4 we 

4 The "down" S-w&ve intensity stays high and nearly constant from m„ w 750 MeV to 
the /o(980) while the "up" solution exhibits a resonant-like maximum under the p(770) . 
5 Hereafter "relatively narrow" means "with a width close to r p = 150 MeV". 
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compare our predictions to the 7r°7r° data and use both sets of data to improve 
our knowledge of the nir S'-wave . The results are summarized in Sect. 5. 

2 Description of the tttt data, first comparison 
and relative normalization 

There were several experiments yielding 7r 7r° data from reaction (3). The 
highest statistics was provided by the GAMS experiment which collected 1.5 x 
10 6 events |3f at 38 GeV/c and 6.4 x 10 5 events 0] at 100 GeV/c. However 
these high-energy experiments were mainly dealing with the high m nn region. 
The NICE experiment ]15| at 25 GeV/c, the KEK E135 experiment [H| and 



some, even older, experiments were of inferior statistics. In this paper we use 
the results of the E852 collaboration [Tj| obtained at 18.3 GeV/c i.e. very close 
to the CCM momentum of 17.2 GeV/c. Additional (and essential) advantage 
of these results is their availability on the World Wide Web where explicit num- 
bers are given. These are intensities of partial waves in Am,™ = 40 MeV bins 
for four intervals of the four-momentum transfer t between -0.01 GeV/c 2 and 
-1.5 GeV/c 2 . 

For a comparison with our results we have merged two \t\ bins i.e. (0.01 4- 
0.10) GeV 2 /c 2 and (0.1040.20) GeV 2 /c 2 into a single one i.e. (0.0140.20) GeV 2 /c 2 . 
This was done by adding the values and their errors linearly and in quadrature, 
respectively. This binning is exactly equal to that in the original CCM paper 
H while in papers I and II we used Am wn = 20 MeV bins and \t\ = (0.005 4 
0.200)GeV 2 /c 2 . 

It should be noted that both experiments have certain problems at low m n7T . 
The spark chamber CM experiment |I| had low and badly controlled acceptance 
for two close tracks; thus it was decided not to publish data below 600 MeV. 
The E852 data (TT2| start right from the threshold but they are somewhat dis- 



turbed by some K® — > 7r°7r° decays which escaped rejection. The upper m nn 
limit of 1800 MeV (1600 MeV) for \t\ > 0.01 GeV 2 /c 2 (\t\ > 0.005 GeV 2 /c 2 ) 
in the CCM data is of lesser importance for our study of ambiguities present 
mainly below the /o(980) . However a warning should be added concerning the 
tc°tt° data at m W7T = (1400 4 1600) MeV. An inspection of Fig. 5D and Fig. 



7D of ref. [12[ shows a substantial G (J = 4) wave below 1600 MeV, which 
seems to be very unlikely at such a low mass. Consequently the 7r°7r° S'-wave 
is probably underestimated in this mass region. 

Unfortunately the E852 collaboration has provided only the acceptance- 
corrected numbers of events and not the absolute cross section which is available 
for the CCM data. Thus we have determined the normalization to our data us- 
ing the strongest Do wave, generally well fixed by the £q moment and saturated 
by the / 2 (1270). The fit in the range m nw = (1000 4 1600) MeV yields the 
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normalization factor = (8.58 ± 0.16) x 10 6 /i5/event in 20 MeV bins. This 
normalization factor was then applied to other waves. The results are shown in 
Fig. 1. 

An inspection of Fig.l shows: 

- a good agreement in the shape of the Dq wave; this shows that this wave can 
be used for normalization, 

- a systematic difference in the Du wave (unnatural spin parity exchange) with 
virtually all E852 points below the 7r + 7r~ data, 

- a similar but much less significant effect in the D^ wave (natural spin parity 
exchange). 

Now let us concentrate on the Du wave. We attribute the difference to 
the assumptions described in the introduction, which were used in the PWA for 
7r°7T but not for the 7r + 7r~ data. This effect was already observed in the CCM 
paper (see Fig. 5 of ref. [§]), when the intensity of the Du wave obtained using 
the data from the polarized target experiment was compared to that from the 
analysis || based on the standard assumptions. The assumptions led to an un- 
derestimation of the Du wave and the overestimation of the least-fixed S- wave, 
both the total intensity and that of the D - wave being fixed by the relevant mo- 
ments. Assuming the same mechanism we have tentatively calculated the excess 
of the 7r + 7r~ Du intensity over the 7t°7t° one and subtracted it from the S-wave 
intensity in the latter. This leads to a (17 ±4)% reduction of the S'-wave in the 
7r°7r° plot and the normalization factor F' s = (7. 13 ±0.35) x 10~ 6 /ib/20 MeV per 
one event. Including the above underestimation of the D^ wave by about 10 % 
yields F'J = (6.44 ±0.52) x 10~ 6 jj,b/20 MeV per one event. Although this whole 
procedure looks like a step in the wrong direction (generally the 7r°7r° data are 
already below the 7t + tt~ ones) we shall see later that a good agreement in one of 
the S-wave solutions can be reached after properly handling the isospin 1 = 2 
component of the S'-wave amplitudes. Since we do not know whether F s or F s 
is the proper value we use the intermediate value F s = (6.70±0.78) x 10~ 6 /ib/20 
MeV. The error range covers both values including their error bars. The value 
F s was used to normalize the E852 data and Fig. 2 shows a comparison of the 
cross sections found in paper I for two solutions: "up-flat" and "down-flat". 
The errors of the 7r°7r° intensities were obtained by adding the normalization 
and statistical errors in quadrature. It should be noted that we have repeated 
the calculation of the normalization factor for various mass ranges as well as 
for the CCM binning (Am„ = 40 MeV and \t\ > 0.01 GeV 2 /c 2 ); the results 
change very little. 

The correctness of the above normalization is supported by the following 
reasoning. Since the tttt production mechanism in the D-wave is very similar to 
that in the S-wave we can estimate the isospin 2 contribution to the D-wave 
amplitude in the 7r + 7r~ and 7r 7r° intensities dominated above 1 GeV by the 
isospin resonance ^2(1270). It was found in Ref. || that in the 7T7T effective 
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1: Intensity of D-waves for the tt°tt° data (full circles) and the 
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Fig. 2: Intensity of S"-wave for the 7r 7r data (full circles) and the 7r + 7r data 
(open circles). Solid line denotes the 1 = 2 contribution 2 x |0.5C 7r t/T 2 | 2 , see 

Eqs. (|27|) and (0). 
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mass range between 1.25 and 1.50 GeV the 1 = 2 D-wave phase shift is negative 
and its absolute value does not exceed 4 degrees. We have estimated the D 
wave isospin 2 contribution assuming that the .D-wave a\ exchange amplitude 
is smaller than the corresponding one pion exchange amplitude and that the 
dependence of the vertex form factors on t is the same for the tt + h~ and 7r°7r° 
production in both S- and D-waves. It turned out that the Dq wave intensities 
in both reactions should be equal within about 2 %. The assumption about 
the vertex form factors is justified by a very similar shape (to about 5 % below 
\t\ = 0.4 GeV 2 ) of the t- distributions in the / 2 (1270) and p(770) mass ranges 
as shown in Ref. 

Let us also comment on the errors being much smaller for the E852 results. 
This is mainly due to high statistics and to smaller numbers of partial waves 
(e.g. the low mass S-wave saturates the 7r°7r° intensity while in the 7r + 7r~ data 
it competes with much stronger P-wave). On the other hand the E852 errors do 
not include the systematic errors of the standard assumptions inherent in their 
PWA. These assumptions were absent in the CCM analysis. In comparison of 
the both sets of data the large normalization errors are very important. 



3 Calculation of S— wave intensity for 7r 7r chan- 
nel from our analysis of the tt + 7T~ data 

In this chapter we shall discuss relations between the tttt production ampli- 
tudes in various charge combinations. In addition to the reactions (1) and (3) 
we consider the following reaction: 

7T + p — ► 7l + 7T + n (4) 

which is an essential source of information about the tttt phase shifts in the 
scalar isospin I = 2 tttt channel M. 

At high energies of incoming pions and at small momentum transfers to the 
target the above three reactions are dominated by the one pion exchange mech- 
anism. However, as shown in || and in I the pseudovector ai(1260) exchange, 
although smaller than the one pion exchange, is also present and cannot be ne- 
glected. Let us denote by T^tttt — > tttt) the off-shell elastic or charge exchange 
mi amplitudes and by t(ai7t — > tttt) similar amplitudes in which the exchanged 
ai is off-shell. The coupling of the exchanged pion to the target proton and 
neutron is relatively well known. This is not the case for the a± coupling but 
in paper I we have been able to separate the a>\ exchange amplitudes using 
the polarized target data of reaction (2). Now we shall use those amplitudes 
and assuming the isospin symmetries of the T and t amplitudes we shall write 
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down the amplitudes corresponding to the 7r°7r° production. Using the isospin 
symmetry |17]] we can write: 



T(7r+7r - _> 7r+7r - ) = 1 Tg + l Tl + 1 T (5) 

6 2b 

T(n+n- -> ttV) = -1t + ^T 2 , (6) 

T(7T + 7r + -> 7T + 7T + ) = T 2 , (7) 

where Tj (J = 0,1,2) are the mi isospin (I = 0,1,2) amplitudes. They are 
normalized to Argand's form: 

Vj e 2 ^ - 1 

Tjf = 2. ' (8) 
where ^ are inelasticities and <5j are the channel j phase shifts. Similarly one 
can write the din — > 7r7r amplitudes as follows 

t(a 1 + 7T- -> 7T + 7T-) = -t + -tl + ~h, (9) 

o 2 o 

t(0i + 7T- -> TrV) = --t + gt 2 (10) 
t(ai+7T + -> 7T+7T+) = t 2 , (11) 

where tj (j = 0,1,2) are again the isospin amplitudes corresponding to the 
production of the tttt pair by the a\ meson scattered on a pion. Let us notice 
that the isotensor amplitudes T 2 and t 2 enter in different combinations in the 
above equations. After the partial wave projection the isospin 1=1 amplitudes 
will not contribute to the S-wave tttt production so the 7\ and t\ amplitudes 
can be ignored. 

As in paper I we consider two independent S-wave irn transversity am- 
plitudes g and h for the reaction (2) which differ by a proton or neutron spin 
projection on the axis perpendicular to the production plane (g =< n j \T\p f> 
and h =< n | \T\p |>). By g and ho we denote the amplitudes corresponding 
to the reaction @ and by g 2 and h 2 the amplitudes of the reaction (|J). 

The amplitudes of the type g or h for those three reactions can be expressed 
as sums of two terms proportional to the proper combinations of the amplitudes 
Ti and U. Thus according to paper I : 

g = C„U{\t q + \t 2 ) + C a V(h + h 2 ), (12) 
6 b 6 b 
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h = C,U\\t q + ~T 2 ) + C a V*(h + ~i 2 ), (13) 
3 b 3 6 

where the functions C n and C a are proportional to the propagators of the ex- 
changed tt and cii , respectively. The expression for the C n factor is following: 



c = -V 4, e 

where: 

p„- is the incoming 7r~ momentum in the u~p cm. frame, 

is the pion mass, 
s is the square of the total energy in the u~p cm. frame, 
q n is the final pion momentum in the tttt cm. frame, 
g 2 /4:TT = 14.6 is the pion-nucleon coupling constant |] 

and / is the complex correction factor. This phenomenological parameter should 
account for the averaged t— dependence of the pion-nucleon vertex function and 
an additional small phase of the pion propagator. In paper I this parameter, 
close to 1, was introduced to satisfy the requirement that the average inelasticity 
coefficient for the S - wave isoscalar tttt amplitude should be equal to 1 below 
the KK threshold. As in paper I we take the slope parameter a = 3.5 GeV -2 . 
The complex functions U and V which depend on the kinematical variables like 
m n7T and t are defined by Eqs. (16) and (17) of paper I. The exact functional 
form of C a (which can be also obtained from paper I) will not be needed in 
further considerations. 

Using Eqs. © and fllOD we can express the amplitudes go and h for the 
reaction (^) in the form analogous to Eqs. (|T2|) and flT3"|): 



9o = C^Ui-^To + ^T 2 ) + C a V(~to + h 2 ), (15) 

h = C«U*(-~T + 1t 2 ) + C a V*{- l -to + h 2 ). (16) 
Finally we give equations for the amplitudes of the reaction (|J): 

92 = C n UT 2 + C a Vt 2 , (17) 



h 2 = C^U*T 2 + C a Vt 2 . (18) 

6 This value of the pion-nucleon coupling constant has been used by us in paper I. According 
to new re sult of Ref. g 2 /47r = 14.17 ± 0.05 ± 0.19. Since CV depends on the ratio 

y (3 2 /47r)// any change in the value of the pion-nucleon coupling constant can be easily 
incorporated by a proper change of the modulus of /. 
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The S"-wave double differential cross section for the reaction (^|) reads: 

d 2<J i 1 2 i 17 1 2 



dm^dt 



In order to calculate it from (jig ) and ( jig) one needs to know four complex 
amplitudes T , T 2 , to an d t 2 . Even if one exploits the knowledge of the ampli- 
tudes g and h from our analysis of the polarized target data for the reaction 
(H) it is impossible to calculate (jl9|) and find direct relations between g and 
go or between h and /io- This would be simpler if the E852 had made an ex- 
periment on polarized target. Since this was not a case we have to make an 
additional assumption to simplify the amplitudes. We shall neglect the ampli- 
tude t 2 . This assumption was made in past in the approximate determination 
of the isotensor-scalar amplitude T 2 M . Now we do the same believing that t 2 
is much less important that the amplitude t which we do not neglect .[] 
Thus the equations (|T2]), (|T3|), ( [T3|) and (|T^) can be written as follows: 



g = C*u£T + ~T 2 ) + C a V(±t ), 
6 b 6 


(20) 


h = C«U\\t q + \t 2 ) + C a V*(\t Q \ 
3 o 6 


(21) 


go = C W U(~T + ^T 2 ) + C a V(~to), 


(22) 


h = C n U*(~T + l -T 2 ) + C a V*{~t ). 


(23) 



Now it is easy to relate amplitudes for two reactions in question: 

go = \c^UT 2 - g, (24) 

h = \c^U*T 2 - h. (25) 

From the above equations we infer that apart of the sign the amplitudes 
g, h and go, ho, corresponding to the 7r + 7r~ and 7r°7r° production differ mainly 
by the presence of the term which includes the 1 = 2 amplitude and enters in 
different proportions in (|5]) and (|j). This fact has an important consequence on 
the rn-x-K dependence of the cross section (|l9f) . 

In practical calculations we need amplitude averages over some ranges of 
the momentum transfer squared t. In paper I we have calculated the am- 
plitudes g and h averaged over —t between 0.005 and 0.2 GeV 2 /c 2 , so over 
At! = (0.2 - 0.005) GeV 2 /c 2 = 0.195 GeV 2 /c 2 . Let us suppose that we need 

7 The standard approach was to ignore both amplitudes to and ^2- In this analysis we 
ignore only t% which is a correction to to- 
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the amplitudes g Q and ho averaged over some other range of —t namely At (for 
example we shall take At=(0.2-0.01) GeV 2 /c 2 = 0.19 GeV 2 /c 2 ). Then one ob- 
vious renormalization factor of the amplitudes is n = (At/Ati) 1 / 2 . This factor 
will be used for the parts of g and h corresponding to the a\ exchange where 
we do not expect a rapid variation with t (for details see paper I). The part 
corresponding to the one pion exchange (proportional to the function CV) is 
more sensitive to the t-variation so we introduce another factor R equal to the 
ratio of the averaged product of C^U over the range At to the averaged product 
of C n U over the range At±. Then using the definitions of the complex factors 
c\ and C2 (Eqs. (32) and (33) of paper I ) one can derive the following relations: 



g Q -=~C 7r UT 2 + b g g + b h h (26) 



and 

^- 1 



where 



h = -C n U*T 2 + b h *g + b g *h } (27) 



, = -{ Cl R + Cl *)n, (28) 
b h = c 2 (l - R)n (29) 



and C n U and C n U denote the averages of the products C W U and C W U* over 
the At range. Using ( PB"D and ([27|) we obtain the effective mass distribution 
integrated over the momentum transfer squared: 

lM> e/ = |^| 2 + |^| 2 . (30) 



We shall also call it the intensity Iq of the 5*-wave in the reaction (^). 

Thus using the analysis of the CCM 7r + 7r~ data performed in paper I and 
knowing the 1 = 2 component we can calculate the Jo intensity for the reaction 
(3) using the phases of the g and h amplitudes fixed in paper I. Appearance 
of the new data on the 7r°7r° production enables us to relax some assumptions 
about those phases. 

In paper I the phases of the transversity amplitudes g and h have been 
obtained essentially from the assumed resonant p(770), ^(1270) and p 3 (1690) 
phases and the measured differences of the S-wave phase and the above men- 
tioned P, D and F waves. However, in order to account for the P and the D 
wave phase differences in the range of m T7r above 980 MeV, we have introduced 
an empirical additional phase A which was linearly decreasing function of 
reaching the zero value at 1420 MeV (see Eq. (52) of I). Below 980 MeV where 
there was not enough available partial waves the value of A was arbitrarily 
assumed to be equal to 50.37°. Now we can directly calculate the values of 
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the A function independently for each m n7T bin using the new 7r°7r° results for 
Jo in addition to the 7r + 7r~ data. In this way one can adjust the phases of g 
and h to fit simultaneously the tc + tt~ and the 7r°7r° cross sections. Let us note, 
however, that using the 7r°7r° intensity as a single input at a given effective tttt 
mass we cannot calculate two independent phases of g and h. Therefore below 
1460 MeV we follow the parametrization of paper I for the phases d g and of 
the amplitudes g and h: 

f tff - 0? + e p{m for 600 MeV < m nn < 980 MeV, 

9 \ tff - df + 0/ 2 (mo) + A for 980 MeV < m nn < 1460 MeV, 1 j 

9 / #h ~ < + ^(770) + A for 600 MeV < m W7r < 980 MeV, 

1 $h - K + %(i270) + A for 980 MeV < m nn < 1460 MeV. 1 ' 

In the above equations & p (77o) and #/ 2 (i27o) are phases of the resonant amplitudes 
defined in paper I and the remaining differences of the $ values represent the 
measured quantities in the partial wave analysis of the CCM group. The nu- 
merical values of A have been calculated to reproduce the experimental values 
of Iq. In subsequent fits we assume that the errors of phases $ g and $h are 
solely given by the experimental errors of the differences between phases of the 
S and P or D waves determined in the partial wave analysis of the 7r + 7r~ data. 

Linear interpolation of the 40 MeV bin data of the E852 group (including 
errors) has been applied in order to obtain the intermediate values corresponding 
to the 20 MeV m nn bins of the CERN-Cracow-Munich group. Following the 
arguments given in Section 2 the effective mass 1460 MeV is the maximum value 
at which we have used the E852 data. 

4 A joint analysis of the 7r°7r° and 7r + 7r data 

Fig. 3 shows the comparison of the 7r°7r° intensity Iq from eq. (|30|) to the 
E852 data using the normalization factor F s determined in Sect. 2. Comparing 
Fig. 2 and Fig. 3 we immediately see an essential role of the 1 = 2 component. 
Now it is the "down-flat" solution which reasonably agrees with the E852 data 
apart from the region above 1460 MeV, for which we suspect an underestimation 
of the >S-wave in the E852 data (see Sect. 2). On the other hand the "up-flat" 
points are too low, especially in the region of m n7T = (870 -j- 970) MeV. We have 
fitted the E852 data to both 7r + 7r _ solutions with the normalization factor f s 
as a free parameter. This has been done in two mass ranges: 

- up to 980 MeV where the solutions are most different, 

- up to 1460 MeV to avoid the uncertain mass region in the E852 data (see Sect. 
2). 
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Table 1 



Table 1: Results of the fits to the data 



Solution 


Mass range (MeV) 


X 2 /NDF 


/, (units lO^Vb) 


"down-flat" 


600 -T- 980 


8.4/18 


7.11 ±0.42 


"down-flat" 


600 1460 


16.6/42 


6.82 ±0.37 


"up-flat" 


600 -T- 980 


39.2/18 


4.33 ±0.33 


"up-flat" 


600 1460 


45.9/42 


4.16 ±0.29 



The results are shown in Table 1. While the fit is very good for the "down- 
flat" solution, the one for the "up-flat" version shows the strong disagreement. 
In the low mass range there is about 3cr deviation^ with the normalization 
coefficient f a being even more significantly (3.4(7 ) different from the F s = (6.44± 
0.52) x 10~ 6 fib/20MeV which is the lowest possible normalization factor (see 
Sect. 2). In the whole mass range the fit is good but demands the normalization 
factor f s which is 3.8a below the corresponding F s . If we had not corrected for 
the influence of the standard assumptions and used the normalization factor 
then the 7r°7r° points would go up by approximately 20%. This could be still 
(although barely) reconciled with the "down-flat" solution but the deviation 
from the "up-flat" solution would be even more significant. 

Let us stress that the fits to the "down-flat" solution yield the normalization 
factor f a between F s = (7.13 ± 0.35) x lO" 6 ^6/20 MeV and F s " = (6.44 ± 
0.52) x 10~ 6 fib/20 MeV in good consistency with both. In particular the fit in 
the (600 -T- 1460) MeV mass range yields both value and error of f s very close to 
F s = (6.70 ±0.78) x 10" 6 fib/20 MeV from Sect. 2.0 This confirms that we have 
applied proper correction to the Du and S"-waves. Now all the D-waves (where 
the 1 = 2 component is negligible) nicely agree between two experiments while 
the S'-wave for E852 data is entirely consistent with the Io for the "down-flat" 
solution of paper I. 

Before concluding on a final rejection of the "up-flat" solution we have 
critically examined the assumptions introduced in paper I and repeated here. 
The aim of the exercise is an attempt to modify some of the assumptions in such 
a way that the "up-flat" solution is saved and we reproduce consistently both 
7T°7r and 7r + 7r~ distributions. An obvious candidate is our assumption 
concerning the A (see eq. (|3l|) and eq. (32)) - phase difference appearing in the 
S-wave g and h amplitudes. 



8 Hereafter a denotes a standard deviation. 

9 This value of f s means raising the E852 Djj wave exactly to the level of the tt + tt~ data 
and a small increase of the E852 Dm wave - see Fig. 1. 
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In paper I we dealt with only four observables from the CCM data in each 
mass bin. These are: 

^ = \d\ 2 + l^| 2 " intensity of the 7r + 7r~ S'-wave , 

\g\/ 1 h\ - ratio of the moduli of transversity amplitudes of the above wave, 
two phase differences (for g and h) between S-wave and the higher wave which 
is resonating in the region in question. 

These four observables were used to determine four real parameters (two 
complex amplitudes) so we had to make some assumptions on A, which enters 
the expressions for the last two observables. Our ansatz assumed the A for 
the S-wave to follow the behaviour of the higher partial waves (P and D) for 
m n7T > 980 MeV. However for lower mass the D-w&ve is too weak to serve as 
a reference wave and we have arbitrarily assumed A to be constant and equal 
to that at m n7T = 990 MeV. As shown above it is in the low mass region that 
the "up-flat" solution disagrees with the E852 data. Therefore we check: 

- whether one can find the values of A providing an agreement of the E852 data 
with the "up-flat" solution as well as with the "down-flat" solution, 

- whether such set of values leads to reasonable behaviour of main parameters, 
namely the inelasticity r] and the phase shift 5. 

This was done by using the intensity Iq of the 7r°7r° S'-wave from the E852 
data in addition to the observables listed above. A minor difference with respect 
to paper I is using the actual resonance parameters from the last Review of 
Particle Physics edition Also the range parameter for the p(770) meson 
width was chosen equal to 5 GeV -1 for both solutions "down" and "up" (the 
Breit-Wigner parametrization of the resonant amplitudes is written in paper I). 

We have critically reexamined all the solutions of partial wave analysis done 
by the CCM group. In case of more than one solution for a given ("up" or 
" down" ) branch we have consequently used the continuity criterion. Since these 
solutions had overlapping errors, the changes are small. At five points of 
(at 630 MeV for the "down-flat" solution and at 990, 1130, 1250 and 1330 MeV 
for the solutions "down-flat" and "up-flat") new values of the CKM data were 
chosen. However, those differences in data have changed the values of 5 and r\ 
by much less than their errors calculated in paper I. 

We proceed in the following way: first the moduli of the transversity ampli- 
tudes g and h are calculated from the first two observables, then the 7r°7r° data 
are used to determine the values of A in each mass bin. We use the normal- 
ization factor F s = 6.70 x 10~ 6 [ib/20 MeV and calculate the vertex correction 
factor / = \f\e l( ^ (see eq. |14|) by minimizing the (rj — l) 2 difference for 19 points 
of m nn < 980 MeV. Now we have: 

\f\ = 0.814 ± 0.019, <p = (-5.4 ± 6.8)° for the "down-flat" solution compared 
to \ f\ = 0.84, <p = -17.8° in paper I, 

|/| = 0.924 ± 0.027, <p = (-32.3 ± 4.3)° for the "up-flat" solution compared to 
|/| = 0.89, <p = -4.4° in paper I. 
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The resulting values of A are compared in Fig. 4 with the ansatz of paper I. 
The variation of A as a function of the effective mass is remarkably rapid for the 
"up" case between 700 and 900 MeV. The values of A drop suddenly from about 
100° to negative values as low as —40°. This jump is related to a considerable 
difference between the Iq intensities seen in the same mass range in Fig. 3. On 
the other hand for the "down" case the values of A are positive and do not vary 
too much in that range. They are also closer to our educated guess on the A 
made in paper I. At 990 MeV the value of A cannot be determined for both 
solutions since the value of Io is particularly low at this effective mass (see Fig. 
3). The same difficulty is found for two points of m n7T greater than 1420 MeV 
and smaller than 1460 MeV. Also at 610 MeV for the "up-flat" solution the 
value of Iq is too large to satisfy Eqs. fl2"BD, (P7|) and (ftOD, so at all the above 



masses there are no corresponding points in Fig. 4. 

Having new / (the same for each mass bin) and new A (different in each 
bin) we proceed with the calculation of four quantities, above all the phase 
shifts 5 and the inelasticities rj of the isoscalar titt — > titt reaction amplitude, as 
well as the phases and moduli of the a\ exchange amplitude. In each mass bin 
we demand them to describe five experimental quantities i.e. the 7r°7r° S'-wave 
intensity in addition to four 7r + 7r~ observables Q used in paper I. The fitting 
program finds solution and errors in each bin, usually with a very small x 2 . 
The exceptions are 990 MeV, 1430 MeV and 1450 MeV bins for both solutions 
and 610 MeV for the "up" case. Here the x 2 is larger (~ 0.1) but reasonable. 
Above 1460 MeV the fits are not acceptable as it can be expected from Fig. 3. 

The results for the most important parameters i.e. the phase shift 5 and 
the inelasticity rj are shown in Fig. 5 and Fig. 6. Unfortunately an interesting 
region around 1500 MeV is hardly covered by this analysis because of the high 
mass problem in the E852 data. Therefore for larger than 1460 MeV we 
show the results of paper I. 

For the "up-flat" solution the inelasticity 77 significantly deviates from unity 
for m n7r = 980 MeV being too large below 700 MeV and too small above 850 
MeV. This behaviour of 77 is unphysical since there are no known inelastic 
processes which can couple to the irn channel (see paper II). Also the phase 
shift 5 decreases substantially between 750 MeV and 900 MeV. The unphysical 
behaviour of 77 and S around 900 MeV reflects the large deviation of 7r 7r° data 
from the predictions based on the 7r + 7r~ "up-flat" solution - see Fig. 3. Thus 
a change of A or / cannot save the "up-flat" solution once the 7r°7r° data are 
included. This conclusion is independent of the normalization uncertainties. 

On the other hand the "down-flat" solution behaves in a reasonable way. 
Comparing this solution in Fig. 6 with our results from paper I (see Fig. 7a 
and Fig. 8) we see that the additional constraint of the E852 data brought 



10 They are: S'-wave intensity, |<?|/|/i|, fig and ■dh- 
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a noticeable reduction of errors. In turn comparing the present solution with 
the standard one of Ref. [1] we note some differences in spite of an overall 
similarity (see Fig. 7b). They are clearly seen in the low mass range. The 
present solution starts at lower 5, then crosses the old one around 700 MeV and 
stays above it up to the rapid change around 990 MeV. 

We have estimated systematical errors due to our particular choice of phases 
fig and ~&h in Eqs. (31) and (32). If one determines the phase -dh as in Eq. (32) 
with A taken from paper I and calculates d g to reproduce the experimental 
values of Iq then the corresponding results for 5 and rj lie always very well 
within the errors displayed in Fig. 6. Thus we conclude that these systematic 
errors are substantially smaller than the errorrs resulting from the joint fits to 
the 7r + 7r~ and 7r°7r°data. 

Let us finally comment on the role played by the pseudovector (ai exchange) 
transversity amplitudes. These amplitudes are significantly different from zero. 
Their moduli are shown in Fig. 9{£] The ratio of the moduli of the pseudovector 
amplitude to that of pseudoscalar (it exchange) amplitudes fluctuates around 
0.2 4-0.3 except of the mass regions around the KK threshold and around 1500 
MeV where moduli of both amplitudes are comparable. In those two regions 

11 Let us note here that in paper I on Fig. 7a) the moduli of the pseudovector amplitudes 
were incorrectly shown due to replacement of the "down-flat" points by the "up- flat" points. 
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there exist two scalar resonances /o(980) and /o(1500). A closer insight into an 
interplay of the a\ and the 7r exchange amplitudes in the effective mass distri- 
bution reveals particularly small values of the interference terms. This result 
may be the main reason while our more sophisticated analysis has produced the 
solution for the isoscalar S— wave somehow resembling the old one [0J outside 
of the ranges near the /o(980) and /o(1500) resonances. 

In our opinion these results represent now the best knowledge on the tttt 
5*-wave in the mass range m n7T = (600 -j- 1460) MeV. They are based on the 
best and 7r + 7r data and use the weakest assumptions about the irn 

amplitudes. 

5 Summary 

We have compared the ir + 7c^ S-wave production data obtained at 17.2 
GeV/c by the CERN-Cracow- Munich collaboration with the 7r°7r° data recently 
collected by the E852 group. This comparison has been made in two steps. In 
the first step we have used the one pion and a x exchange model developed in 
paper I to predict the S-wave intensity of the n n production at 18.3 GeV/c. 
Here the previous assumptions about the phases of two independent g and h 
helicity amplitudes have been directly taken from paper I. Presumably small 
isotensor part of the a\ exchange amplitude has been neglected. The proper 
normalization of the S*-wave 7r°7r° intensity has been done using the strong 
.Do~waves in both experiments. The result of the comparison of the calculated 
7r°7T S-wave intensity with the 7r 7r data was satisfactory for the "down- 
flat" solution but the 7r°7r° intensity corresponding to the "up-flat" solution 
was significantly and systematically too low in the m nn range between 800 MeV 
and 1000 MeV. 

In the second step we have varied the phase correction A to the g and h 
amplitudes at each m n7T bin trying to fit simultaneously the tt + tt~ and 7r 7r° S- 
wave intensities. Then we have recalculated the tttt scalar-isoscalar phase shifts 
and inelasticities. Again for the "down-flat" solution to the 7r + 7r~ data we have 
obtained a reasonable behaviour of the inelasticity parameter which, as one 
should expect, was close to unity below the KK threshold. However, the m n7T 
variation of inelasticity corresponding to the "up-flat" solution was too strong. 
Also the "up-flat" phase shifts have an unphysical minimum at 900 MeV. This 
is a clear signal that the "up-flat" solution should be rejected. Since the "down- 
flat" phase shifts and inelasticities behave regularly we accept the "down-flat" 
solution, resulting from the joint analysis of the n + ir~ and 7r°7r° data as a unique 
physical one. 
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